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Abstract 

We investigate the renormalization group (RG) flow of SU(3) lattice gauge theory in 
a two coupling space with couplings Pu and /3i2 corresponding to 1 x 1 and 1x2 loops 
respectively. Extensive numerical calculations of the RG flow are made in the fourth 
quadrant of this coupling space, i.e., /3ii > and /3i2 < 0. Swendsen's factor two blocking 
and the Schwinger-Dyson method are used to find an effective action for the blocked gauge 
field. The resulting renormalization group fiow runs quickly towards an attractive stream 
which has an approximate line shape. This is a numerical evidence of a renormalized 
trajectory which locates close to the two coupling space. A model fiow equation which 
incorporates a marginal coupling (asymptotic scaling term), an irrelevant coupling and 
a non-perturbative attraction towards the strong coupling limit reproduces qualitatively 
the observed features. We further examine the scaling properties of an action which is 
closer to the attractive stream than the currently used improved actions. It is found that 
this action shows excellent restoration of rotational symmetry even for coarse lattices with 
a 0.3 fm. 



1 Introduction 

Since Wilson's first numerical renormalization group (RG) analysis of SU(2) gauge theory O, 
there have been many Monte Carlo RG studies of non-perturbative /3-functions (see Ref.[0] 
and references therein). In these analysis indirect information about the (3 function, such as 
A/3, has been obtained [^. Recent progress of lattice techniques [@, |5|, ^] allows us to estimate 
directly the RG flow in multi-coupling space [0. 

We study renormalization effects by means of a blocking transformation which changes 
the lattice cut-off but leaves the long range contents of the system invariant. A new blocked 
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action S' as a function of blocked link variables V^s is constructed from the original action 
S{U) as 



-S'{V) 



J P[7e~^(^)5(y - P{U)), (1) 



where P defines the blocking transformation. The action S' includes the renormalization 
effects induced by blocking. In the space of coupling constants, the blocking transformation 
makes a transition from a point corresponding to 5 to a new point, S' . Repeating the 
blocking transformation, we obtain trajectories in coupling space which define the so called 
renormalization group flow. 

There is a special trajectory, i.e. renormalized trajectory (RT), which starts at the ultra- 
violet fixed point. On the RT, the long range information corresponding to continuum physics 
is preserved. Recently Hasenfratz and Niedermayer have stressed that the action on the RT 
can indeed be considered as a "perfect action" Q. Therefore if we find a RT corresponding 
to a blocking transformation, it provides an action which gives accurate results corresponding 
to the continuum limit. Even if it is an approximate one, it serves as a well-improved action. 
In this sense, a pioneering work has been done by Iwasaki more than ten years ago [^. 
He estimated a RT by matching Wilson loops based on a perturbative approximation, and 
proposed an improved action which we will call below Iwasaki action. 

In this work, we analyze numerically the RG flow in two coupling space, (/3ii,/?i2)) of 
SU(3) lattice gauge theory and clarify the structure of the renormalization group flow. The 
action is restricted to the following form; 



plaq rect 



(2) 



Here Ppiaq and Prect correspond to 1 x 1 and 1x2 loops, respectively. 

The main purpose of the present work is to perform an extensive study of the RT beyond 
a perturbative analysis. It is a non-trivial fact whether the RT locates near the two coupling 
space, i.e. in the (/3ii,/3i2) plane. If no remnant of the RT can be seen in this plane, the two 
coupling space is insufficient to obtain good improved actions. In this sense, a global analysis 
of the RG flow from weak to strong coupling regions is indispensable. Parts of our analysis 



were reported in Refs.[^, 17 1, and in this paper we present much more data which are enough 
to confirm the existence of RT. 

Our analysis is performed in the fourth quadrant of the coupling space. We examine 
renormalization effects induced by Swendsen's factor two blocking on the plaquette action as 
well as on some improved actions. In addition, we try to clarify the global structure of the RG 
flow. An evidence that the RT sits close to the two coupling space is provided by the fact that 
the flow runs quickly towards a narrow attractive stream. Characteristic features of the flow 
in the strong and weak coupling regions are also found. The observed features are reproduced 
by a model flow equation which incorporates a marginal coupling (asymptotic scaling term), 
an irrelevant coupling and a non-perturbative attraction towards the strong coupling limit. 

Based on the flow structure, we examine the scaling properties of several actions defined 
in this two coupling space. Tests are made for the rotational invariance and the scaling of 
\fa jTc- Near the attractive stream, we find good restoration of the rotational invariance. 

This paper is organized as follows. In sect. 2, the basic tools for the analysis are given. 
Sect. 3 is devoted to present simulations and numerical results of the RG flow. Scaling tests 
are described in sect. 4. In sect. 5, a model flow equation which can reproduce the observed 
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RG flow is proposed. Renormalization effects beyond the two coupling space are discussed in 
sect. 6. A summary of results is given in sect. 7. 

2 Blocking transformation and determination of renormaliza- 
tion effects 

Here we describe the basic framework to study the RG flow for SU(3) gauge fields on the 
lattice. First we produce field configurations with an action S which has coupling constants 
(/3ii,/3i2, • • •)) s-iid apply a blocking transformation on these configurations. Next we deter- 
mine an action S' with coupling constants {P'u, Pi2r ' ') which reproduces the transformed 
configurations. In this way, we extract a flow in the coupling constant space from the original 
coupling constants {(3} to the transformed ones {/?'}; this is the RG flow. 

In order to perform the blocking transformation on a lattice, we adopt Swendsen's factor- 
two blocking [|ll|. For a set of SU(3) link variables U^{n), a blocked field is constructed 
as: 

Qt,{n) = U^{n)U^{n + /i) + c ^ Uy{n)U^{n + v)U^{n + i> + fi)U^^{n + i> + 2/i). (3) 

Here c is a parameter to control the weight of the staple- like paths. A convenient nota- 
tion V-fj,{n) = Vil{n — ft) is also used and the sum is taken over negative as well as posi- 
tive direction. Q^{n) is projected onto the blocked SU(3) gauge field Vf^{n) by maximizing 
ReTr(Q^(n)y;(n)). Q 

To determine the effective action S' on blocked configuration V, we use the Schwinger- 
Dyson method Q . This is based on the following identity: for a link V/^ , consider the quantities 

{lmTiiX'Vi,Gl)) = \j'DV Im Tr {\%,Gl) e"^' (4) 

where A** stands for Gell-Mann matrices. Here action is assumed to have the form, 

S' = Y. Eyll-^R^Try.G/] (5) 

I 7 

and CfJ stands for a "staple" coupling to the link I in a loop of type 7. For the present analysis, 
7 corresponds to a plaquette and a rectangle. Eq. (^) should be invariant under the change 
of variable V/^ — > (1 -|- ie\^)ViQ. Setting terms linear in e to be zero, we get the identity, 

j VV[ReTT{{X'fVi,Gf^) + ImTr(A^l^„G^,)ImTr(AVz„G,„)] e^^' = 0, (6) 

where 

Gi = Y.^-^G]. (7) 

7 

Summing over b in the expression (|6|) above, we obtain the Schwinger-Dyson equation, 

^ The projection from Q onto a SU(3) matrix V is not unique. In ref. j^, we employed the polar decompo- 
sition, i.e., Q = VH where H is an Hermite matrix and V = V/det(V') with V = QHQ'^Qf''^- One can easily 
prove that both methods are equivalent for unitary matrices Vs. 
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+Re(Tr(Gr„(G;:,)t)) - ^(Im Tr(l1„GrjIm Tr{Vi,Gl))}. (8) 

Here we have used the identity: J2b=i Tr{X^A)Tr{X^B) = 2Tr AB - fXrATrB. We apply eq. 

to the blocked configurations, and calculate the expectation values (• • •) on both sides. 
Now eq. (^) may be considered as a set of linear equations with /3^'s as unknowns. It is noted 
that we may use other loop operators instead of G". However a minimal choice is to take the 
same G°'s as the ones entering the action. In this case, the number of equations is equal to 
the number of unknown couplings. 

It is noted that the canonical Demon method also works for the present purpose 
This method tunes the effective action so as to reproduce the mean values of the plaquette 
and rectangular loops whereas the Schwinger-Dyson method respects wider loops which are 
combination of staples such as Tv{Gf{G^)^). In case of a limited coupling space, they may 
lead to different actions. Because the systematic errors involved in both methods are not 
known in the present stage, only the results obtained via the Schwinger-Dyson method are 
presented in this work. 

Since we study the KG flow in the two coupling space, we show corresponding Schwinger- 
Dyson equation explicitly. 



[(3[, ) [ A22 J [ {ReTriPji'J)) J ^ ' 

where 

All = i^E.^^Re[(TV(P«P«t)^ _ (TV(P«P«)) - i(Tr(P«)TV(p(;.)t 
A12 = ^E.^^Re[(Tr(p(;.)pg)^)) - (TV(P«P(^;)) - i(TV(P«)TV(p(f - 

+ (Tr(P«P^^;)) - (TV(P«P^J) - i(TV(P«)TY(P^^; - Pp)] 
A21 = ^,E.^^M{Tt{P^I^pL'J^)) - (TV(P£)P«)) - i(TV(p£))TV(P«^ -Pi^)))] 
A22 = ^E.^^Re[(Tr(p(^.)p(f )) - (TV(P£)P£))) - i(TV(p(^.))TV(p£)^ - p£))) 

+ {Tr{P^l^P^J)) - (MP^l^Pp) - 1{T,{pLI^)Tt{P^J - Pp)] 

with 

P^l\n) = V^{n)V^{n + fi)V-f,{n + // + iy)V-^{n + u) 

P£^(n)= V^{n)Vuin + n)Vu{n + n + u)V-^{n + n + 2u)V-^{n + 2u)V-uin + iy) n-^-s 
Pjluin) = V^{n)[Vy{n + ^l)V.^,{n + /u + v)V^^,{n + u)V^y{n - /i + u)V^{n - fi) 
+Vi,{n + iJL)Vy{n + 2/x)y_^(n + 2^ + v)V-^,{n + + v)V-y{n + v)] . 



3 Simulation and numerical results of the RG flow 

Using the techniques in the previous section, we study the coupling flow induced by the factor 
two blocking. We set the blocking parameter c = 0.5. 

Series of simulations on lattices of size 8^ and 16^ are performed. The region surveyed in 
the present work is the fourth quadrant of the (/3ii,/3i2) plane which covers most improved 
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fii2 




7.00 0.0 
7.00 -0.35 
11.00 -0.9981 
13.20 -1.493 
(t)13.20 -1.493 


13.188(22) -1.6564(94) 
8.148(29) -1.0431(98) 
15.189(73) -2.329(19) 
15.445(50) -2.559(15) 
15.519(37) -2.583(13) 



Table 1: Examples of renormalized couplings fS'n and f3i2 obtained through the Schwinger- 
Dyson method. A hundred configurations are used at each point. For (f ) , measurements are 
performed every 1000 pseudo heat-bath steps, while for the rest of the data we measure every 
100 steps. 



(in (3i2 


Pll Pl2 


lattice size 


6.20 0.0 


9.515(56) -1.156(23) 
9.547(14) -1.1671(53) 


8'* deconfined 
16^ confined 


9.8496 -0.8937 


12.368(61) -1.883(16) 
12.280(14) -1.8688(41) 


8^ deconfined 
16^ confined 


9.0 -1.03332 


7.211(38) -1.046(17) 
7.2086(86) -1.0456(30) 


8^ deconfined 
16^ confined 



Table 2: Renormalized couplings (3[i and (^[2 on the deconfined and confined lattices. These 
are calculated with 50 configurations, and measurements are performed every 100 steps. 

actions presently known. We use the pseudo heat-bath method to generate the gauge fields. 
The blocking transformation is carried out at more than 30 points. At each point, about 
100 configurations separated by every 100 sweeps are used to determine the renormalized 
couplings (3[i and f3[2 through the Schwinger-Dyson method. Examples in the determination 
of the coupling are shown in Table |^ and The errors are given by the Jack-knife method 
and they are relatively small even at deconfined points. Analyses are made also for data 
samples in which configurations are separated by 1000 sweeps. Those data agree with that of 
the standard samples within error bars. An example is shown in the last two lines of Table |l|. 

We check also the effect of the finite temperature phase transition on the determination 
of the coupling constants. Since the blocking transformation induces renormalization effects 
corresponding to a change of lattice cutoff, the resulting coupling shifts are insensitive to the 
phases. We examine this point by comparing results belonging to the confined phase on the 
16^ lattice while remaining deconfined on the 8^ lattice at several values of f3 as shown in Table 
^. The finite temperature phase transition on the 8^ and 16^ lattices for the plaquette action 
takes place at approximately f3u = 5.9 and 6.3 which roughly correspond to a^/a = 0.27 and 
a^/a = 0.14 respectively. As seen in Table |2|, both results agree within errors. 

Our results are summarized in Fig. |^, in which the coupling shift resulting after the 
blocking is indicated by arrows. 

As shown in the figure, several characteristics of the flow are seen. If we start from the 
plaquette action {P12 = line), renormalization results in a negative Pu as expected by a 
perturbative analysis. At Pn = 6 ~ 8, the renormalization effect is very strong making Pn 
twice larger and P12 negative. The resulting points are far below the line of the tree-level 
Symanzik action. On the other hand, in the strong coupling region below Pn < 5, the 
effect of the renormalization is to reduce Pu. Therefore, the plaquette action suffers large 
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Figure 1: Renormalization group flow of SU(3) lattice gauge theory in the two couphng space 
iPii^Pu)- There is an attractive stream to which the arrows converge. Dotted and dashed 
hnes correspond to ( tree level ) Symanzik action and Iwasaki action respectively. The long 
dashed line corresponds to the "DBW2" action introduced in section 4. On the DBW2 line, 
the black dot is the point reached by twice blocking from the plaquette action at (3 = 6.3 on 
a 32^ X 64 lattice 



renormalization effects and it is far from the renormalized trajectory. 

Points starting on the line defined by the tree-level Symanzik action, Pu/Pu = —0.05, 
are RG transformed onto ones with (3i2 far more negative. Although renormalization effects 
are reduced, the trend is the same as that of the plaquette action. This means that tree-level 
Symanzik action is still far from blocking invariant. This indicates that perturbative O(a^) 
improvement is insufficient at least for the currently used range of lattice spacings. 

Now we look at flows starting from points on the line corresponding to Iwasaki action, i.e., 
Pi2 = — 0.09073/3ii . We see that renormalization occurs approximately along the line up to an 
intermediate point (/Jn « 7.3, /?i2 = — 0.09073/3ii). At larger Pn, however, flows depart from 
the line defined by this action. Thus the present blocking transformation renormalizes Iwasaki 
action further and induces more negative values of P12 above the intermediate coupling region. 

Let us turn to the global structure of the flow. Blockings are made starting from points 
with /3i2//3ii = —0.1 ~ —0.15. For those points, the renormalization effect is relatively small 
and converges to a narrow stream. This trend is manifest at strong coupling. As a whole, there 
is an attractive stream which the flow approaches quickly. Furthermore, once the flow reaches 
the stream, it runs along it. Therefore actions on the attractive stream are approximately 
blocking invariant apart from a normalization. The shape of the attractive stream is clearly 
recognized as a parabolic curve in the strong coupling region while at points far from the 
origin, it is less obvious with the present data. This is a remarkable indication that the 
renormalized trajectory locates close to the (/3ii,/3i2) coupling space. This is an encouraging 
result for finding a good improved action in this two coupling space. 

A closer look at the attractive flow allows to extract more information. If we start from 
the Wilson action at Pn > 6.0, the flrst blocking leads to larger renormalization effects as Pu 
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increases and the resulting flow vectors have the same direction pointing towards increasing (3ii 
further until the flow has reached the main attractive stream. This feature seems consistent 
with a flow induced by an irrelevant coupling. In the strong coupling region the behavior 
is quite different, under the blocking the coupling moves deeper into strong coupling, with 
the main stream following a parabolic behavior, as already indicated above. In section 5, we 
will try to reproduce these features by a model equation including an irrelevant coupling, an 
asymptotic scaling term and a driving term derived from the area law behavior of Wilson 
loops. 



4 Tests of improvement for actions near the attractive stream 

The purpose of this section is to examine the degree of improvement for actions lying near the 
attractive stream. As shown in Fig. ||, we expect that the irrelevant coupling is considerably 
reduced and actions will be dominated by the marginal coupling term. Here we analyze an 
action suggested by double blocking from Wilson action at /? = 6.3 on a 32^ x 64 lattice where 
{f3n,f3i2) = (7.986(13), -0.09169(41)) has been obtained 0. This action is called DBW2 
(doubly blocked from Wilson action in two coupling space ) and has a value of /3i2//3ii equal 
to -0.1148 (long dashed line in Fig. |^). 

There are many ways to test improvement. Among them, we will examine rotational 
invariance of the heavy quark potential and independence of Tc/ ^/a on the lattice spacing. 
Let us define a measure of violation of rotational symmetry as 



where V{R) is the static quark potential and 5V{R) is its error. Von{R) is a fitting function 
to only on-axis data. X^off means summations over off-axis data. This quantity is measured 
on the configurations generated by the DBW2 action. For comparison, plaquette-, tree- level 
Symanzik and Iwasaki improved actions are also examined. A 12^ x 24 lattice is used for 
this purpose and simulations are made for lattice spacings ranging from a = 0.15 to ~ 0.4fm. 
The statistics is a hundred configurations for each data point and the error is given by the 
Jack-knife method. 

The results are summarized in Fig. ^ 

In the figure, the horizontal axis is the lattice spacing squared so as to see the expected 
O(a^) violation. As seen in the figure, both Iwasaki and DBW2 actions show excellent restora- 
tion of rotational symmetry even at a ~ 0.3fm while clear a? violations are seen for plaquette 
and tree-level Symanzik actions. 

A second check is the scaling of Td \fo. The critical temperature is defined as 

T, = XjNta, , a, = a(/3e) (13) 

where Nt is the temporal lattice size. In order to obtain the critical coupling f3c , Nt = 3, 4, 6 
lattices are used and the Polyakov loop susceptibility is measured. Here, histogram method 
is utilized to determine the peak of the susceptibility, [p^] Then, [3c at infinite volume limit is 
obtained by finite size scaling of 12^ x 4 and 16^ x 4 lattices. The resulting /3c are given in 
Table |. 

The string tension is measured on 12^ x 24 and 18'^ x 36 lattices at the values of the 
coupling determined by the Tc analysis. It is extracted from the static quark potential using 
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Figure 2: Rotational symmetry violation 6v versus for various improved actions. 



Nt 


Ns 


+ 8/3i2)c at y = oo 


3 


10, 12, 14 


0.75696(98) 


4 


12, 16 


0.82430(95) 


6 


18 


0.9636(25) 



Table 3: Phase transition points for the DBW2 action at finite temperature. Results are 
obtained by extrapolating to the infinite volume limit. 

the Ansatz ; 

V{R) =A + - + (jR 
Jx 

From the extracted values of a at Pd^t = 3), (3c{Nt = 4) and (3c{Nt 
following values of Tc/ -y/a: 

Tc/^/a = 0.6340(60) at Nt = 4 

0.6301(65) at Nt = 6 (15) 

In Fig. ^ we show the obtained results together with other data from plaquette (Wilson) , tree- 
level Symanzik, tadpole improved Symanzik and Iwasaki actions taken from Refs []l9| , pc| ]. 

No appreciable dependence on the lattice spacing is seen in the ratio Tc/\/a for all the 
cases including the DBW2 action. 

The results of both tests is an indication that the DBW2 action is a well-improved action 
as suggested by the non-perturbative RG flow analysis. 



Pu + 8/3i2 


A 


a 


•> 

aa^ 


0.82430 
0.9636 


0.550(17) 
0.5791(96) 


-0.255(23) 
-0.357(22) 


0.1555(28) 
0.06996(99) 



Table 4: Parameters of the heavy quark potential for the DBW2 action. These simulations 
are carried out on 12'^ x 24 and 18'^ x 36 lattices while the values of the coupling are those 
determined by Tc analysis (Table |3|). Statistics is 130 configurations per data at Pu + 8P12 = 
0.82430 and 80 configurations per data at Pu + 8P12 = 0.9636. 
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6) we obtain the 
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Figure 3: Scaling behavior of Tc/^ct. We compare the DBW2 action with results from other 
actions in two-coupling space |19, 20 1. 



As a by-product of the Tc analysis, contours of constant lattice spacing are obtained 
through the relation a = l/NtTc- Fig. Q shows a compilation of the phase transition points 
in the two coupling space for the plaquette, tree-level Symanzik, Iwasaki and DBW2 actions 
on A^t = 3 ~ 12 lattices. 
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Figure 4: Contours of constant aTc in the /9i2)-plane. Symbols represent the phase 
transition points for each Nt at the thermodynamic limit. Dotted lines connect the points 
with the same aTc{= ^/Nt) and they suggest contours of aTc = const.. 



5 A model analysis for the flow 

Results of the RG flow in the section 3 provide basic information on relation between couplings 
and lattice spacing. In order to understand the driving force of the flow, we examine a toy 
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model which incorporates a perturbative scahng term and an irrelevant coupling at small 
lattice spacing and an attraction originated from area law in strong coupling region. 

We consider two dimensional model beta function which has a marginal coupling and an 
irrelevant one at small lattice spacing as 

^ = Ap-B{Pn)T^ (16) 
ama |n|^ 

where B{x) is the perturbative beta- function 

Bix) = 126o + 72- + • • • (17) 

X 

with bo = 33/(487r^) and 6i = (102/121)6q. A is a constant 2x2 matrix which has a zero 
eigenvalue with an eigenvector w and a negative eigenvalues A with a eigenvector v . n is 
constant vector orthogonal to v , i.e., (n-v) = . By this condition, resultant solution exibits 
correct scaling property for {n ■ fj). 
Solution of eq. ([l6[) is 

^ = {ca^ + d{a))v + '^^^w (18) 
[n ■ w) 



where /3p(aA) is the asymptotic scaling solution satisfying 

-B{^) (19) 



d(5p 



din a 

and 



d{a) = i?:^ /' di r'-'Bi(3p{CaA)) . (20) 
(n -w) Jo 

It is noted that {n ■ (3) = Pp{aA). 

Attractive driving force in the strong coupling region comes in from the string tension. 
The lowest order of strong coupling calculation for N x M Wilson loops is (suppose that NM 
is an even integer) 

< WiN X M) >= (/3n/18)^^^ + (/3ii/18)^^^"2(^^^/^8)piVM ^ 

(/3n/18)^^^-n/3i2/18)2p2^^ + • • +(/3i2/18)^^^/2p^i^/2 (21) 

where are the tiling weights for filling the area of the loop by {NM — 2k) [1 x 1] and 

/c [1 X 2] tiles. Then, area law for the Wilson loop leads 

/3ii/18 = exp(-aV), /3i2/18 = C exp(-2aV) . (22) 

In a differential form, we have 

Driving force in a two dimensional toy model beta function is assumed as sum of those in 
the weak and strong-coupling region,i.e.. 



10 



Here, we keep the next order in a in the non-perturbative term as free parameters, Ci and ^2- 
Parameters included are 6,9', XXi and (2- v and iB are specified as iT = {cos 9', sin^') and 
m = {cos 9, sm9). Here the matrix A is A = AiT (8> vjt^/{v ■ wr) with wt = i(72W. n is given 
by n = (l,cot(-6l')). 

By this model, we calculate lattice spacings for different actions, RG flow trajectories and 
contours of constant lattice spacings. Comparisons with the present data and those reported 
in references |]28| , [29(1 are made. Giving initial conditions for fin , (5i2 and a on the (3i2 = line 
(plaquette action) 23, 22] , those quantities are calculated and shown in Fig.s |^, and 
Parameters are chosen so as to match the recently reported data of a^/a for the Iwasaki 
action [28| 



-0.156 , 9' 



-0.205 , A 



-0.5 , Ci = 0.1 C2 



0.0 



(25) 



As shown in the figures, lattice spacings in the — /3i2 plane and the RG flow trajectories 
are reasonably well reproduced. As for the flow, rapid approach to the attractive stream is 
driven by the irrelevant coupling. At intermediate lattice spacings, both the non-perturbative 
and asymptotic scaling terms drive the trajectories. Finally, in the strong coupling region, 
the non-perturbative term dominates and the parabolic behavior of the trajectories sets on. 
We also note that the model describes contours of constant lattice spacing fairly well as seen 
in Fig.|6[ 

Through the study, an understanding for basic driving mechanism of the RG flow in 
— (312 plane by scaling term and an irrelevant coupling at small lattice spacing and an 
attraction originated from area law in strong coupling region is suggested. 





Figure 5: (a) Lattice spacing versus fin. Crosses are input data from the plaquette action 
given by ||2l|, Filled symbols are the data of tree-level Symanzik action (squares) 

and DBW2 action (diamonds) ||2^. The corresponding open symbols are the results pre- 
dicted by eq. (p^. (b) Flow trajectories calculated by eq. (|2^). The trajectories start 
at (/3ii,/3i2) = (5. 55, 0.0), (5. 70, 0.0), (6. 20, 0.0) and (6.80,0.0). Symbols on the trajectories 
indicate the points corresponding to the factor two blocking where the lattice spacing is 
2"'ao(ri = 1, 2, 3, • • •) with initial value oq. 



fa =420MeV is assumed for the data in reference 123 
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Figure 6: Contours of constant lattice spacing. Filled circles show the points with a^/a = 
0.4099 , 0.2702 and 0.1619 as computed by eq. (M). The data for {aTc)''^ = 4 (open squares) 



and 6 (open circles) obtained by the scaling analysis (Fig. ^ are also indicated. 

6 Examination of renormalization effects outside of the two 
coupling space 

In this work, the working space is limited to two coupling space (/?ii,/3i2) • We have seen an 
evidence that the RT locates close to that plane. Since this is a highly non trivial fact, the 
truncation effect which comes from the limitation of the coupling space to the two-dimensional 
plane should be examined. There are several indirect indications that truncation effects may 
be important. As seen in Fig. |l], the shape of the attractor is less clear above ~ 12 . This 
suggests that the distance to the RT grows in this region. 

Another point is that the resulting couplings {P'u, (3'i2) obtained through the truncation 
into two coupling space after blocking give larger lattice spacing (roughly 40-60 %) than that 
of the original blocked lattice. One possible reason is the choice of a minimal set of loops 
in the Schwinger-Dyson equation. But this effect might also be due to the truncation of the 
space of couplings. See ref. p^ . 

In order to examine renormalization effects outside to the two coupling space, we perform 
a limited analysis in a three coupling space, (/3ii, /3i2 5 Awist), with a twisted loop included in 



the action |2(:] 



•^twist — Awist Ln.At^y^rll - 5^eTr(P/^''igJ] (26) 

where 

PZlt = U^{n)Uy{n + fi)Ur{n + (l + u)U-^{n + ft + u + -f)C/_^(n + fi + f)[/_,(n + f). (27) 

We will study Awist in addtion to (/3ii,/3i2) because in Ref. |17| more general actions were 
studied and it was found that the value of the twist coupling is larger that that of the chair 
type loop. Moreover the action with VFixi, Wix2 and VFtwist is often used in the reteratures. 
1^, We perform the blocking transformation starting from points near the attractor in 
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Figure 7: RG flow of SU(3) lattice gauge theory in three coupHng space /3i2, Awist)- 
Blocking is applied for the points near the attractive stream on the /3i2)-plane at /3twist = 
0. The top figure shows the projection of the flow onto the Awist) -plane while the bottom 
figure is the projection onto the (/3ii,/3i2) -plane. 

the (/3ii,/3i2) plane, i.e. the /3twist = sector. The results are shown in Fig. 0. As shown 
in the figure, no sizable value of /3twist is generated for points on the attr actor. On the other 
hand, if we start from the points ~ 10,/3i2 ~ —1.0 , a negative /3twist is induced. This 
indicates that the attractor in three coupling space sits below the /3twist = sector in this 
region. Although more extensive studies are necessary to clarify the exact situation, it seems 
that the attractor in two coupling space gives a good starting point for designing improved 
actions. 

7 Summary 

We investigate the renormalization group (RG) flow of SU(3) lattice gauge theory in a two 
coupling space, (/9ii,/3i2)- An extensive numerical calculation of the RG flow on the lattice is 
made. Swendsen's blocking followed by an effective action search using the Schwinger-Dyson 
method is adopted to flnd renormalization effects. In our previous analysis, we adopted the 
Demon method. Although both methods give similar results, the Schwinger-Dyson algorithm 
is more echonomical because in case of the Demon formula we need an additional Monte Carlo 
simulation for each conflguration. Moreover the Schwinger-Dyson method uses larger loops 
to determine the couplings. 
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Analyses arc performed in the fourth quadrant of the coupUng space and reveal the pres- 
ence of an attractive stream. Trajectories are first attracted towards this stream and after- 
wards they move towards the origin along it. The stream converges to a parabolic curve in the 
strong coupling region. These features indicate that the RT locates close to the two coupling 
space and the attractive stream traces the RT. 

A model flow equation which consists of asymptotic scaling, an irrelevant coupling and a 
non-perturbative force corresponding to the area law can reproduce the observed features. 

In this paper we have compared several actions in the two coupling space by measuring the 
restoration of the rotational symmetry and the scaling of ^fajTc- In the regions of o ~ 0.3, 
Iwasaki and DBW2 actions, which are near to RT , are superior to tree-level Symanzik and 
plaquette actions. Although " improveness" of all actions are indistinguishable at smaller 
lattice spacing within the present statistics, this indicate that nonperturbative study of RG 
flows is valuable to design an improved action at intermediate lattice spacing. It is highly 
desirable to pursue the same analysis for fermion actions. 

The effect of truncation to a space with only two couplings is partly examined and it is 
found that renormalization effects outside the two coupling space are small and the attractive 
stream in this space gives hence a good starting point for improvement. 
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